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Abstract 

We consider the spin-exchange in ultracold collisions of hydrogen (H) and antihydrogen 
(H) atoms. The cross sections for transitions between various hyperfine states are calcu- 
lated. We show that the process of spin exchange in H — H collisions is basically driven 
by the strong force between proton and antiproton, that spin-exchange cross-sections are 
proportional to the difference of singlet and triplet strong force S-wave scattering lengths, 
while the strong force effects are enhanced about 10 9 times due to long range atomic forces. 
We show that any difference in hyperfine structure of H and H ( CPT violation ) manifests 
itself in radical change of spin-exchange cross-sections energy behavior. 
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1 Introduction 



The successful production of low energy antihydrogen atoms [TJ [2] and current experiments 
on their trapping [3, 4 J have stimulated interest in the interaction of antiatoms with ordinary 
matter. The simplest example of matter- antimatter interaction is given by collisions between 
antihydrogen and hydrogen atoms [5j [6j [7l [9l [T(J [Til HI] ■ Collisions between antihydrogen 
and helium atoms have also been studied [T3] . These investigations revealed that cold collisions 
between atoms and antiatoms do not necesserily result in annihilation, but may lead to elas- 
tic scattering or rearrangement reactions. Interestingly, it was shown that the cross sections 
for processes that do not result in annihilation (e.g. elastic scattering) are sensitive to the 
strong interaction between particles and antiparticles [Ml [13l [TSl [12l [TBI HZ] ■ This presents an 
interesting perspective for investigating the details of strong interactions between matter and 
antimatter in low energy collisions of atoms and antiatoms. 

In the present work we investigate spin exchange in hydrogen-antihydrogen collisions. The 
oncoming atoms are assumed to be asymptotically in the well defined hyperfine states which 
can be changed as the result of collisions. This can happen with or without changing the 
internal energy of the individual atoms so that some of these collisions appear elastic, apart 
from spin-flip. We find that in the case of H — H scattering the main source of such spin- 
exchange transitions is the strong force. This is in contrast to spin-exchange in atom-atom 
collisions (see [TSl 112] and the references therein), where the strong- force effects are usually 
negligible. At the same time the leptonic contribution to spin-exchange in H — H collisions is 
found negligibly small, in spite of the fact that the magnetic moments of the leptons (electron 
(e) and positron (e)) are on the order of 10 3 times larger than those for the hadrons (proton (p) 
and antiproton p)). We will show that long range van der Waals interactions between atoms and 
antiatoms plays role of magnifying glass that exposes the strong force effects on the molecular 
scale. Another unique feature of the spin-exchange cross sections is their sensitivity to the 
equality of the hyperfine energy levels in H and H , ensured by the CPT invariance. Violation 
of that invariance results in a radical change of certain cross sections for spin-exchange. 

The paper is organized as follows. In section[2]we present the formalism used in calculating 
spin-exchange cross sections. Section [3] discusses the main results. The physical sense of these 
results is analysed in terms of a transparent zero-range potential model of strong forces. Section 
[4] contains the conclusions and the Appendix gives the details of certain derivations. 

2 Formalism 

We assume the following four-body Hamiltonian for the HH system: 



Here H nr is the nonrelativistic spin-independent Hamiltonian, which includes the kinetic and 
Coulomb interaction energy for all particles. The subscript a labels the pairs of interacting 
particles: {pp), (pe), {pe), (pe), (pe) or (ee). W a (s a ,r a ) is the spin-dependent pair interaction, 
arising from low order relativistic corrections to the nonrelativisic Hamiltonian; W a (s a , r a ) acts 
on the spin s Q and the relative coordinates r a of the particle pair a. 

The forthcoming analysis can be greatly simplified by taking into account the characteristic 
spatial and energy scales of the interactions involved in H — H scattering. As pointed out in 
[201 [211 [22l [5j [23l [6l [7l [24J, [25] , the rearrangement transitions to Protonium (Pn) and Positronium 
(Ps) occur mainly at internuclear distances smaller than the rearrangement radius R r ~ 1 a.u.. 
For internuclear distances above R r the adiabatic approximation, leading to the separation of 
leptonic and hadronic motion, is well justified. This allows the description of ultralow energy 
H — H scattering within a one-channel model based on a nonrelativistic Hamiltonian [8 , 23 , 24 : 




(1) 



a 



H nr ~ Tkin + Vopt + Vad(R)- 



(2) 
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The above Hamiltonian describes the interatomic motion in a fixed adiabatic leptonic state, 
where T^ in is the operator of kinetic energy of internuclear motion, V a d{R) is the adiabatic 
potential in a given leptonic state, and V opt is the effective complex optical potential that 
accounts for the inelastic rearrangement and vanishes for interatomic distances R > R r . 

In the following we will discuss the ultralow energy collisions of H and H in the ground 
leptonic state. We will use the adiabatic interaction potential V a d(R) obtained in 26J. It is 
worth noting that since the Pauli principle is not forbidding an electron and a positron to occupy 
identical spin-space states, the adiabatic potential V ac i{R) is independent of the leptonic spin 
state. This is in contrast to the case of HH, where due to the Pauli principle the adiabatic 
potential is different in singlet and triplet electronic states. This difference is important for our 
considerations. In fact, the spin-exchange transitions in H — H collisions can occur only when 
the (relativistic) spin-dependent terms W a (s a ,r a ) are taken into account. At the same time in 
the HH case the main contribution to such transition rates comes from the above mentioned 
effective spin dependence of the nonrelativistic adiabatic potential, while the explicitly spin- 
dependent interaction appears as a perturbation. We will show in the following that the spin- 
exchange reactions in H — H collisions are particularly sensitive to the spin dependent strong 
force between proton and antiproton. 

The general form of the spin-dependent interaction is: 

W a = W° a {r a ) + W^(r a )(L a S a ) + W* a (r a )(r a S a f + W* b (r a )S 2 a . (3) 

In the above formula r a , L a , S a are relative coordinate, angular momentum and spin operators 
for the particle pair a. 

The particular form of the lepton-lepton and lepton-hadron pair interactions can be obtained 
within the QED [27] . In the case of (pe) and (pe) they are responsible for the hyperfinc splitting 
of the ground state of hydrogen and antihydrogen. In the present work we will approximate 
the non-relativistic pair interactions for (pe) and (pe) by the following effective terms |18j : 

WW(ep) = a H F(s e ,eip,p) = -7p(F 2 - 3/2). (4) 

Here cihf = 2.157 10~ 7 a.u. is the hyperfine constant, s ej g is the electron (positron) spin, i Pt p 
is the proton (antiproton) spin and F = s + i is the total spin of the (pe) or (pe) pair. The 
above interaction correctly reproduces the ground state hyperfine splitting of separated H and 
H atoms. 

The influence of the (pe) or (pe) spin-dependent terms on the spin-exchange reaction rate 
is expected to be small due to the Coulomb repulsion that prevents the involved particles from 
approaching each other. We will neglect such terms in our study. 

In the case of (ee) the spin-dependent interaction is responsible for the fine structure and 
annihilation in Positronium [27]. Due to the CP invariance this interaction conserves spin in 
Positronium. In our treatment we will use the following simplified model of the (ee) spin- 
dependent term: 

W (es) = a°d(r Ps )P Q + a^irp^Pt. (5) 
Here a and a 1 are the singlet and triplet interaction constants: 

a° = -4.3 • 1(T 4 - i4.9 • 1(T 6 a.u., (6) 
a 1 = 3.3 • 1CT 4 - i4.3 • 1(T 9 a.u., (7) 

Po = \S Ps = 0,0)(S Ps = 0,0| and P\ = J2m\Sps = l,M){S Ps = 1,M| are the projection 
operators on the subspace of singlet or triplet states of Positronium. The suggested form of 
interaction conserves spin and correctly reproduces the energy shift and annihilation lifetime of 
para- and ortho-Positronium. The spin-orbit coupling terms are neglected; we will show below 
that such terms are not important for the spin-exchange transition rates in the H — H system. 

The leading contribution to the spin-dependent component of the (pp) interaction arises from 
the strong force. Strong interaction is CP invariant and conserves the spin Sp n of Protonium. 
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Furthermore, this interaction is attractive and localized to distances within R < R s with R s ~ 1 
fm (1.88 10~ 5 a.u.). It also includes an absorptive component that is responsible for coupling 
to the (pp) annihilation channels. 

There are several models of low energy pp interaction( see [28] and references therein) that 
differ in the detailed behavior of spin-dependent potentials. These models, however, give close 
values for the pp strong force singlet and triplet scattering lengths a^ 1 . Since the range of 
strong forces is much shorter compared to that of atomic interaction, it is possible to take into 
account the strong force effect through an appropriate modification of the boundary condition 
for the scattering wave function. Such a procedure is based on a knowledge of the strong force 
scattering lengths only [TTJ [16] . Indeed, at distances such that R s <C R <C R r the interaction 
is dominated by the Coulomb pp potential. The wave function of internuclear motion at such 
distances behaves as a superposition of regular F (pR) and irregular G (pR) Coulomb solutions, 
with momentum p — > |16j . The coefficient that mixes these two solutions, the i^-matrix 
element k, is determined by the strong force scattering length according to [SU] HH] : 

fc = tan(5 s °' 1 ) = -2 7 rMa°' 1 ! (8) 

where M stands for the reduced mass of pp. We will use the following values of the 5-state 
strong force scattering length [301 HH] : 

a° sc = (1.07-il.45)-lQ- 5 a.u., (9) 
a\ c = (1.68-tl.Q6)-lQ- 5 a.u. (10) 

The corresponding nuclear fc-matrix elements become 

tan(£° c ) = -0.06 + iQ.08, (11) 
tan(<5] c ) = -O.l + zO.06. (12) 

In terms of state vectors \f(R), Sp n ), R s < R <§C R r , this result can be formulated as follows: 

l/OR), S Pn ) ~ \F Q ( P R),S Pn ) + k\G (pR),S Pn ) (13) 
K = -2nM (a a sc n Q + a^i) • (14) 

Here the 7To,i a re projection operators on the subspaces of singlet and triplet Protonium states. 
In this way the effect of strong forces could be taken into account by introducing the if-matrix 
(I14p . The latter determines the boundary condition that can be imposed on the wave function 
in the range R s < R <§; R r . The effect of the strong force induced boundary condition (fl4j) is 
equivalent to the action of the zero range potential [31j : 

V. = g (<&7? + a s x c 5fi) ^( R )(^^)- (15) 

To proceed with the calculation of spin-exchange rates within the effective adiabatic model 
one needs to find the contribution of all spin-dependent terms ([3]) to the effective internuclear 
interaction. We will accommodate these terms to the first order of perturbation theory. To 
do so we average the spin-dependent interactions that depend on leptonic coordinates over 
the unperturbed leptonic wave function ^ e S (R, r e , r g ) (calculated without the account of spin- 
dependent interactions) . Taking into account our assumptions about the spin-dependent terms 
and their radial dependence, such averaging will involve only the (ee) interaction: 

V eS (R)= J \* et e(R,r e ,r s )\ 2 W eg dr e dr s (16) 

V ee (R) = A(R) (a°P Q + a'PA , (17) 



where A(R) = |\I/ ej g(i?, r e = r g )| 2 is the leptonic coalescence probability. We use the values of 
the coalescence probability calculated in ref. |llj . The corresponding spin-dependent leptonic 
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Figure 1: Real part of the spin-dependent leptonic potentials for the singlet (dashed line) and 
triplet (solid line) states. 

potentials for the singlet and triplet states are shown in Fig. 1. We note that the spin-dependent 
interactions ([3]) also modify the final states of Pn and Ps in the rearrangement channels, so 
that the latter may contribute to the effective optical potential V opt - Such a contribution only 
arises at the level of second order perturbation theory and will not be taken into account here. 

We may now formulate the scattering problem for spin-exchange transitions. The spin state 
of the colliding H — H system can be described in different representations, in terms of 16 spin 
basis vectors. Since the spin-dependent terms for (pp) and (ee) interactions conserve spin of 
the corresponding pairs, they will be diagonal in the basis set with given (pp) and (ee) spins. 
Such basis set will be denoted \S( PtP ),M( PtP ),S( e ^),M( e ^). We will hereafter refer to this basis 
as the ^-representation. 

The asymptotic Hamiltonian, which describes the separated H and H atoms and includes 
the hyperfine terms W( pe ) and W( pS -), is diagonal in the basis with given (pe) and (pe) spins Fh 
and Fg, respectively. Such basis set will be denoted \Fjj, Mh, F h , M h ) and will be called the 
^-representation. The scattering state vector \tp) is a superposition of basis spin states: 

\<p)= f{R) FBl M B ,F Bl M B \F H ,M H ,F H ,M H ) (18) 

F h ,Mh,F h> Mr 

The i?-dependent expansion coefficients play the role of hadronic wave functions for the HH 
system in a given spin state of four particles. The asymptotic form of functions f(R) at large R 
gives in a standard way the ^-matrix in the given spin-state representation. The desired wave 
functions are found from the solution to the coupled equation system: 

Y,(^'\H -E\ f i)F(R)^ v = 0, V = 1,2,. ..,16 (19) 

Here F(R)^> ^ is a 16 x 16 solution matrix whose indices label the outgoing (//) and incoming 
(/i) channels, respectively; and r\ labels the sets of spin quantum numbers Fh, Mh, Fh, Mh- 

While solving the coupled channel equations we utilize the fact that for interatomic distances 
R ~ 1 a.u. and below the characteristic values of the local adiabatic potential are much greater 
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than the hyperfinc splitting. Thus even in the zero-energy limit the hyperfine terms may be 
neglected in the solution (TT5)) at short and intermediate internuclear distances. At such distances 
the equation system can be decoupled in the \S( p Mu,p\, >S(e,e)> M^ e ^) basis: 

Id 2 - \ 

+ V opt + V ad (R)+A(R)a s ^ -E ) F{R)^ = 0, (20) 



2M dR 2 



where /i denotes a set of quantum numbers SVg) , ^i(p,p) , S( e ,e) i ^(e,e) ■ Such an equation system 
should be supplied with boundary conditions, imposed at arbitrary R within the range R s <C 
Ro <^C R r , and specified by the strong force if -matrix (fl4|) : 



F(Ro)w = <W [FoipRo) - 2Tra S J pp) G (pRo)) (21) 
F'(Ro)w = <W ( F o(pRo) - 2ira S s ^ ) G' Q (pR )) . (22) 

Instead of solving the equation system (|20| with the boundary conditions (|21II22[) we use the 
more practical method developed in 16J. It amounts to introducing specially fitted boundary 
conditions at the rearrangement radius R = R r : 

(F'(R r )F- 1 (R r )) iM!iM ,=S fl ^p(R r )cot(S + S^ pf) ) (23) 



where p(R r ) = y/2MV a d{R r ) is a classical local momentum given at the internuclear separation 
Rr. The phase shift S takes into account the effect of the optical and adiabatic potentials 
accumulated at distances below R r . The account of strong forces is achieved by adding to the 
phase shift 6 the spin dependent strong force phase shift S S c PF} ■ The additivity of phases is 
explained by the fact that contribution of the optical potential to the phase shift is vanishing 
in the region R < R s , characteristic for strong forces. The phase shift 6 is calculated using the 
model nonlocal optical potential [16) and is given by: 

6 = 0.74 + i0.32. (24) 

The advantage of the above method is that once the phase shift S is calculated, the effect 
of rearrangement transitions is incorporated through the boundary condition (|23p obviating 
cumbersome calculations with nonlocal potentials. The phase shift approach gives a useful 
universal tool for the treatment and analysis of rearrangement transitions. 

Once the (decoupled) equation system is solved in the S'-representation, the solutions (valid 
in the limited range below R ~ 1 a.u.) should be transformed back to the asymptotically 
correct i^-represcntation. This is done by means of a unitary transformation |19j : 



F(R) = U + F(R)U 
U = (S {P:P) , M (PtPh S {eiS) , M (e , s) \F H , M H , F n , M R ) 



(25) 
(26) 



Thus obtained solution F(R) should be matched at some point Rq ~ 1 a.u. with a solution 
<f>, in which now the hyperfine energy splitting is taken into account: 



E 



i d 2 

' 2M dR 2 



+ V ad {R) <W + [U + V (eS) U}^. p + . 



F'iR^F-^Ro) 



= (27) 



[^(Ro^-'iRo)},. (28) 



Here Q^> , P is a diagonal (16 x 16) hyperfine energy matrix, which gives the threshold energies 
in channels with different combinations of total spin Fh and Ffj. Its diagonal values are: 
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Figure 2: The spin-exchange cross sections (r d& ^ ad (solid line) and <r da ^ cd (dashed line). 



C if (F H = 0,F a = 0) 
Qw = a HF \ 1 if (F H = l,Fg = 0) or (F H = 0,F n = 1) (29) 
{ 2 if (F H = l,Fg = l) 

From the conservation of the projection of total spin (M H g — Mr + Mg) it follows that 
the $-matrix has a block structure, consisting of a 6 x 6 matrix ( M = 0) , two 4x4 matrixes 
{M H g = 1 and M H g = -1) and two lxl matrixes (M H g = 2 and M H g = -2). 

The such obtained solution matrix has the following asymptotic form as R — > oo: 

-» exp(-ipR) - exp^pii)^ 1 / 2 ^ 1 / 2 . (30) 

Here p is a diagonal matrix whose elements are p MiM = y/2M(E — <2 M , M ); exp(=Fipi?) is a 
diagonal matrix with elements [exp(=Fipi?)] M:At = exp(=Fzp MiM i?); and £ is the desired S-matrix. 

3 Results and discussion 

We consider the H — H collisions in the absence of an external magnetic field (B = 0). We 
follow the usual notation for the hyperfine states, i.e. the H state \Fjj = 0, Mjj = 0) is denoted 
as a, the state \F H — 1,M H = —1) is denoted b, the state \F H = l,M H = 0) is denoted as c 
and the state \F H = 1, M H = 1) is denoted d. The corresponding spin-states of H are denoted 
a, b, c and d. We present three types of typical behavior of spin exchange cross sections as a 
function of collision energy. The spin-exchange cross sections tr da ^ a j and c da ^ c j are shown on 
Fig. 2; the cross section (7 c( i^da i s shown on Fig. 3. 

According to the CPT invariance, the energies of the da and ad states are equal and so the 
reaction da — > ad takes place even in the limit of vanishing collision energy, E — > 0. In that 
limit the cross section tends to a constant value, Oda^ad{E) 0.04 a.u. 2 . 

The reaction da — > cd can only occur above the threshold energy Q da ^ c d = a HF- The 
threshold behavior of the associated cross section is clearly seen on Fig. 2. Again, according to 
the CPT invariance, the cross sections for conjugated reactions are the same, giving <J da ^ c d = 
G ad^dc- The inverse reaction cd — > da is exothermic (in our case the energies of the final 
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Figure 3: The spin-exchange cross section o c ^ dn . 



states for exothermic reactions are smaller than the energies of the initial states by curf or 
2a hf) so that the cross section varies as 1/v, where v is the incident channel velocity. Such 
behavior is shown in Fig. 3. Hence, for collisional energies E < auF the cross sections of 
exothermic reactions are the largest among spin-exchange cross sections. Let us mention that 
a hypothetical violation of CPT invariance would manifest itself in a radical change of the 
spin-exchange cross sections. In particular, the energy difference between the hyperfine levels 
of H and H (induced by the CPT violation) would result in the appearance of new reaction 
thresholds. Consequently, the cross sections of the reactions da ^ ad would follow the above 
mentioned (v and 1/v) threshold pattern, instead of tending to the established constant value 
in the limit E — > 0. 

For increasing collision energies, i.e. when E ^> uhf, all spin-exchange cross sections 
tend to the same limit. In such a case the hyperfine splitting can be neglected even in the 
asymptotic states. Therefore the S'-matrix is diagonal in .^-representation. The S-matrix in the 
F-representation (needed in our calculation of cross sections) is obtained by the transformation 
specified in eq. (|26p . One can show (using the explicit form of the transformation matrix U [19j ) 
that the S-matrix elements in the F-representation for the reactions da — > ad and da — > cd, are 
connected to the matrix elements in the S-representation in the following way: 

_ ((o,i|S|o,i) + (i,o|S|i,o)) (MjglM) , . 

^da^ad — ^ ^ 1 

_ ((0,1|S|0,1)-(1,0|S|1,0)) " 

D da-*cd — 1 ■ K aZ ) 

Here (i,j\S\i , j) is the S'-matrix element in the S'-representation, with S p p — i and S e g — j 
= 0,1). It follows from |3T1 l32|) that the equality of cross sections (Jda^ad = a da->cd 
requires: 

(1,0|S|1,0) = (1,1|S|1,1). (33) 

The above equality means that the elastic scattering in the state \S p p — l,S e g = 0) is 
the same as in the state \S p p = l,S eg = 1). This can hold if the contribution of the (ee) 
spin-dependent potential (fT7|) is negligible. A calculation made without the inclusion of the 
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potential (jTTJ) fully confirmed this assertion. The scattering lengths aio and a n in the states 
\Spp = 1) &ee — 0) an d \Spp — 1, S eS — 1) turned out to be: 

arc = 5.665 -VI. 216 a.u. (34) 
an = 5.677 -«2.216a.u. (35) 

while the scattering length ai, calculated without inclusion of the potential (|lTf) turned out to 
be: 

a\ = 5.666 — i2.216 a.u.. 

The small difference between these scattering lengths proves that the role of spin-dependent 
ee interaction is negligible. To understand this we note that the ratio \V eS (R)/V a d(R)\ never 
exceeds 10 -3 and reaches a maximum in the range 3 a.u. < R < 5 a.u.. At these distances the 
WKB approximation is still justified. The semiclassical phase accumulated at these distances 
by the wave function (without account of the V eg (R)) is: 



'WKB 



R-2 



y/2M\(Vad{r)\dr. 



The change of that phase due to the account of spin-dependent leptonic potential V e g(R) is : 



Aowkb = / — -dr < 10 — - — 

J Rl VWUW\ 2 

which explains the small contribution of the leptonic spin-dependent interaction to the phase 
of the wave function, and therefore to the S- matrix [16] . Thus we come to the conclusion that 
the rate of spin-exchange in H — H collisions is determined by the pp strong interaction. 

The short-range character of strong forces enables the factorization of nuclear and atomic 
interactions and simplifies the study of the energy dependence of the spin-exchange cross sec- 
tions. At short distances, where spin transitions take place, the energy difference between the 
final states can be totally neglected. Using the explicit form of the transition matrix U (|2"6"| 
and the zero-range pseudo-potential (|15p one can get the transition amplitude between different 
spin states in the distorted wave approximation: 

f a0 = ^k_^£ J *0*(R)*(R)(JLiJ)*0 {K)d3R = <^k^ m a{0) (36) 

Here ^(R-) is the wave function corresponding the initial (final) spin-state a (/3), calculated 
without taking strong forces into account. (We note that the exact expression for the transition 
amplitude can be obtained from ([36|) upon replacing a(R-) by the wave function calculated 
with the account of strong forces.) The only contribution to the spin-exchange amplitude arises 
from the 5*-wave, since the higher partial wave functions tend to zero at the origin. The spin 
exchange cross section is then given by the following expression: 

*af, = ^ |o ' e ~ oSe|2 |^(0)^(0)| 2 , (37) 

Pa ^ 

where p a jj is the channel momentum. It is convenient to perform a further analysis in terms 
of the Jost function fj. Following the definition of the S'-wave Jost function [32] we may for 
our case write 

RjjyPa) 

where cfP a (R) is the regular S- wave radial solution of the Schrodinger equation (neglecting strong 
forces) in spin-state a. The behavior of the regular solution at the origin is 

0°(je-o) = i2 
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so that 



Equation (f37|) now takes the form: 



fj(Pa)' 



P0 

&aB = K 

Pa 



(39) 



(40) 



2fj{p )f.j{p a ) 

Near the threshold the momentum dependence of the Jost-function is known to be [32] : 

fj(pa) = fj(0) exp(-ip a a) (41) 

where a is the atom-antiatom scattering length (without the account of strong forces). Numer- 
ical calculations give the following values for a and l/|/j(0)| 2 : 



l/|/j(0)| 



5.2 -il. 
29067. 



a.u. 



(42) 
(43) 



We can exploit the smallness of the energy thresholds (determined by the hyperfine constant 
oihf) in different spin states, and get for the near-threshold behavior of the spin-exchange cross 
section: 



cr a a = 7T— exp(2 Im a(p a + pp)) 

Pa 



2/3(0) 



(44) 



It follows from the above formula that, in the case p a = pp, the cross section in the zero- 
energy limit tends to the constant value: 



ctq 



2/3(0) 



0.04 a.u. 



(45) 



while in the presence of the energy threshold the behavior of the cross section is determined by 
the imaginary part Im a of the atom-antiatom scattering length: 



Ta/3 



Pb_ 

Pa 



exp(2Ima(p Q +pp))<ro- 



(46) 



Comparison of the above expression, obtained in the distorted wave approximation, with the 
results of numerical calculations, shows that the formula (|46|) is accurate to within a few percent 
for collision energies below 10 -5 a.u. 

The behavior of the cross section (|4"6"|) near the threshold has a clear physical meaning. Spin- 
exchange can only occur when the nuclei approach distances characteristic for strong forces, i.e. 
1 fm. Due to inelastic rearrangement transitions into the Pn + Ps channels the flux of H — H at 
such distances is damped by the factor exp(2Ima(p a +p@))- The damping rate is given by the 
imaginary part of scattering length for H — B collisions. We note that the imaginary part of 
the scattering length can be extracted from the energy dependence of the spin-exchange cross 
sections [33] . 

In the case of exothermic reactions (e.g. cd — > da or bd — > ad) there is an energy excess 
in the final channel: AEp = naHF, n = 1,2. It follows from (f46|) that the cross sections 
for such reactions behave like d n /v in the limit E — > 0. Taking into account that as E — > 
pp — * s/2MnotHF, we obtain for d n : 



V2M 



nctHF 



which gives: 



M 

di = 
d 2 = 



exp(2 lmay2MnaHF)ao, 



8.0 x 10~ 7 a.u. 2 

1.1 x 10~ 6 a.u. 2 



(47) 



(48) 
(49) 
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The ratio of the near-threshold spin-exchange reaction rates bd — > aa and cd — > da is given 
d 2 /di = V2exp (2(V2 - l)lmay/2Ma H F) = 1.375. (50) 



As one can see, this ratio is determined by the product Imay/2AIauF- 

The cross sections of the inverse reactions (da — > cd or ad — ► bd) show a characteristic 

threshold behavior. In the vicinity of the threshold they behave like pd n / \/2Mna.HF- 

It was already mentioned that for E 3> otn f all the spin-exchange cross sections tend to the 

same limit. Within the zero-range potential model of strong forces, this limit is given by: 

a = exp(4Imap)<7o. (51) 

A remarkable feature of the obtained results is the factorization of the spin-dependent strong 
force effect and the effect of long-range atomic interaction (|33]l . While the first appears through 
the difference between singlet and triplet strong force scattering lengths, the second is given by 
the Jost function for the interatomic motion of the H — H pair. This latter factor (|l//j(0)| 4 ) 
strongly enhances (approximately 10 9 times!) the effect of strong forces, due to focusing of the 
flux oip — p towards the center by the molecular potential. The large value of the enhancement 
factor is determined by the large ratio between atomic and nuclear scales. 

Under certain conditions, the enhancement factor |l//,/(0)| 4 could become even larger. This 
additional enhancement is due to the existence of the near-threshold quasi- molecular H H states 
|16j . Indeed, such quasi-bound states manifest themselves as zeros of the Jost function fj(p) in 
the complex momentum plane. The presence of these zeros in the threshold vicinity leads to a 
small value for fj(0). To study the question in more detail we will use the analytical expression 
for the Jost function based on the WKB approximation for interatomic motion at distances 1 
a.u. < R < 5 a.u. (see Appendix): 



1 2irMa 



\fY KB (0)\ 2 sm 2 (7r/8 + n + S) 



(52) 



In the above equation, ao is the characteristic scale for van der Waals interaction between H 
and H [IS]: 

a = \/2MC 6 ^ 3//4 ) ~ 4.99 a.u. (53) 
V 2^2r(5/4) V ' 

Cq = 6.5 is the van der Waals constant for HH interaction, Q, is the semiclassical phase, 
accumulated in the distance range R > R r 

poo 

fl= y / 2MV ad (R)dR = 19.38 (54) 



and S = 0.74 + i0.32 is the complex phase shift (I24[) accumulated at distances less than R r . The 
value of expression (|52|) . obtained within the WKB approximation, turns out to be 



r w*bm = 24830 ( 55 ) 



and is comparatively close to the exact value (|43|) . 

Following the method developed in [16] , we will vary the total phase of the wave function 
ip = tt/8 + 8 + O. The physical sense of such procedure is the following. It is known that 
interatomic interactions can be modified by means of external magnetic or electric fields (see 
[5HI55] and references therein). Without specifying the physical method of such modification we 
may describe the resulting effect as a variation of the accumulated phase of the wave function. 
We will be interested in such modification of interatomic interaction that results in the increase 
of the enhancement factor (|4"5]) . This factor ([4"5]) would become infinite (the Jost function 
fj(0) = 0) if there is a quasi- molecular state exactly at the threshold. In such a case ip = irk, 
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Figure 4: The zero-energy limit of the cd — > dc cross section as a function of the phase variation 
that expresses a modification of the interatomic interaction due to the presence of an external 
field. 



k = 0, 1, 2.... As long as there is an imaginary component of the total phase tp (which represents 
absorption of the H — H flux into rearrangement channels) such a condition can not be satisfied 
exactly; however the enhancement factor will show maxima corresponding to values of the total 
phase satisfying Re tp — irk. As was shown in [16] this condition corresponds to the appearance 
of decaying quasi-bound state in the spectrum of the HH system close to the dissociation 
threshold. 

In Fig. |4] we present the zero-energy limit of the cross section for the reaction cd — > dc as 
a function of phase variation. For certain values of the total phase tp the spin-exchange cross 
section increases more than 50 times in comparison with its "unbiased" value (tp = 20.51 + i3.2) 
obtained in the absence of any external perturbation. Strictly speaking, in such a resonant case 
the distorted wave approximation (|36| is not justified and one should substitute the approximate 
value of ^(0) by the wave function calculated with the account of strong forces. Nevertheless, 
this result indicates that the modification of the interatomic H—H interaction by external fields 
could in principle strongly increase the spin-exchange cross sections, exposing the effect of strong 
forces. Such additional enhancement might also occur in collisions with other atoms, due to 
the sensitivity of the cross sections to the reduced mass of the colliding nuclei [16] . 

4 Conclusions 

We found that the spin-exchange transitions in H — H collisions are caused mainly by the strong 
force. This mechanism of spin-exchange is a unique feature of atom-antiatom interaction, with 
no counterpart in atom-atom collisions. The corresponding transition amplitudes are propor- 
tional to the difference between singlet and triplet strong force scattering lengths. Meanwhile, 
the cross sections acquire molecular dimensions due to the large enhancement of the strong force 
effect; this enhancement originates from the presence of long range, atom-antiatom interaction. 
We have shown that the enhancement factor is on the order of 10 9 . 

The energy behavior of the spin-exchange cross sections turns to be sensitive to the difference 
between hyperfine energy levels in H and H (which might occur due to CPT violation). If 
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such violation were to occure the threshold behavior of the spin-exchange cross sections would 
be radically different from the CPT invariant case. We present a universal description of 
the influence of external fields on the H — H scattering, which suggests that field-induced 
modification of the interatomic interaction can additionally enhance (up to 50 times) the spin- 
exchange cross sections, exposing the strong force effects. 
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5 Appendix 

The present Appendix outlines the derivation of the analytical expression for the Jost function 
fj(Q) for a zero collision energy of the H—H system. We start with the expression for the S- 
wave scattering length for H—H , obtained in 16J under the following assumptions. At distances 
R > Rr ~ 1 a.u. the H — H interaction is well described by the one-channel, adiabatic potential 
Vad(R)- At distances R> Rh ~ 5 a.u. this potential can be accurately reproduced by the van 
der Waals interaction —Cq/R 6 , with C$ = 6.5. There exists a domain Rj- < R < R^ <C R v dw, 
with R v dw — \/2MC& ~ 10.5 a.u., where the WKB approximation is well justified. 

Matching the WKB wave function to the exactly known asymptotic zero energy solution of 
the Schrodinger equation with the —Ce/R 6 potential, one gets the following expression for the 
scattering length in H — H collisions: 

a = a Q (l + cot(^ + Q + sfj (56) 

where ao is 

r(3/4) , . 

ao ^^27frwi)" 4 - 99a - u - (57) 

and £1 = 19.38 is the WKB phase (|54[) . accumulated in the domain i? r < R < oo. 

To calculate the Jost function, we determine the variation of the scattering length (fB"6"|) with 
the small variation of phase A<5: 



sin' 



Aa = —a — ~—. -. (58) 



In the case of the zero-range potential 

Vzrp = ^a sc 6(R)(-^R) (59) 
the variation of phase is given by © : 

AS = -2nMa sc . (60) 

The corresponding variation of the scattering amplitude, induced by the zero-range potential 
(|59|) . can be obtained within the distorted wave approximation: 

M f 

Af = —— J *o(R)*Vz fl p*o(R)d 3 i? - -a sc |* (0)| 2 (61) 



12 



where ^o(R-) is the wave function calculated without the account of the potential (|59[) . We 
now use eq. (|39p to express the wave function in terms of the Jost function and obtain 



a/ = "W- (62) 

The S-wave scattering length is connected to the zero-energy scattering amplitude by the rela- 
tion a = — /(0). Inserting this relation into eq. (I62|) and comparing (|62|) with |58|) . we finally 
obtain the following expression for l/|/j(0)| 2 : 

1/l/ ^ (0)|2 = ■ STZ^Y ^ 
sin (tt/8 + il + o) 
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